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' Abstract 

■ We point out a synergy between T-conjugated oscillation channels in the determi- 

CN ■ nation of the neutrino mass hierarchy with oscillation experiments with relatively short 

! , baselines (L < 700 km), where the matter effect is small. If information from all four 

J> \ oscillation channels Un — > v e , t>n — > P e , z/ e — > z/,j and z/ e — » is available, a matter 

effect of few percent suffices to break the sign-degeneracy and allows to determine the 
£NJ ■ neutrino mass hierarchy. The effect is discussed by analytical considerations of the 

relevant oscillation probabilities, and illustrated with numerical simulations of realistic 
, experimental setups. Possible configurations where this method could be applied are 

the combination of a super beam experiment with a beta beam or a neutrino factory, 
' or a (low energy) neutrino factory using a detector with muon and electron charge 

identification. 

D 

X< ■ 

> 

X ' 1 Introduction 

The determination of the neutrino mass hierarchy is one of the most interesting open issues 
in neutrino physics. Present data allow for the two possibilities normal hierarchy (NH) and 
inverted hierarchy (IH), which are conventionally parametrized by the sign of the difference 
of the mass-squares of the first and third neutrino mass eigenstates: Am^ > for NH and 
Amjj < for IH. Identifying which of these two possibilities is realized in nature is of great 
importance for our understanding of the neutrino mass mechanism, the relation of neutrinos 
to the charged fermions, and the problem of flavour in general. 

Nevertheless, the determination of the sign of Am^ turns out to be experimentally chal- 
lenging. The most promising way seems to be to explore the matter effect in neutrino 
oscillations [1-3]. This can be done with long-baseline experiments [4], atmospheric neu- 
trinos [5], or supernova neutrinos [6]. Alternative methods to determine the hierarchy, not 
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based on the matter effect, have been proposed [7,8], but they turn out to be extremely 
challenging, if not impossible in practice. 

The usual strategy to determine the mass hierarchy in long-baseline experiments is to 
consider configurations where the matter effect is strong, with baselines as long as possible, 
ideally several 1000 km (see, e.g., ref. [9] for a recent work), and neutrino energies close to 
(or above) the resonance energy [3], which for typical densities of earth matter and lAmU ~ 
2.5 x 10~ 3 eV 2 is around 10 GeV. The basic idea is to observe the effect of the resonance, which 
occurs for neutrinos in case of NH and anti-neutrinos for IH. Hence, the discrimination of the 
two hierarchies is based on the difference of the matter effect in CP-conjugated oscillation 
channels, e.g., v e — > and v e — > oscillations. The ultimate setup to explore this effect is 
certainly a neutrino factory with one or two baselines of several thousand km, see ref. [10] 
for a recent study. Other options are intense super beam experiments with a detector at a 
baseline L > 1000 km [11,12], or very long-baseline high gamma beta beam experiments [13]. 

In this work we will consider a different strategy, and discuss the possibility to determine 
the mass hierarchy with experiments in the regime of small matter effect. We will consider 
experiments with baselines in the range 100 km < L < 800 km and in the energy range 
of 300 MeV < E < few GeV. In this regime the size of the matter effect is typically a 
few percent. We point out that if information from all four CP and T-conjugate oscillation 
channels — > u e , — > P e , v e — > i/^, and v e — > is available such small effects can be 
explored efficiently in order to determine the neutrino mass hierarchy. The basic observation 
is that the so-called sign-degenerate solution [14,15], which prevents the determination of 
the mass hierarchy, moves in opposite directions in the plane of #i 3 and 5qp for T-conjugate 
channels. This effect has been observed in ref. [16] in a simulation of CERN based beta 
beam and super beam experiments. Here we discuss the underlying physics and point out 
the general principle without specializing to a specific experimental configuration. This 
method could be applied for example for the combination of generic beta beam (V e — 
and super beam (V^— >V e ) experiments. Furthermore, a low energy neutrino factory [17] 
with a detector (or detectors) capable of muon and electron charge identification would 
offer a place to apply this method. If charge identification for electrons is not possible, the 
neutrino factory (providing the v e — information) could be combined with a ^—^v e 
super beam experiment [18]. 

The outline of the remainder of the paper is as follows. In sec. [5] we discuss the pro- 
posed hierarchy determination by considerations of the relevant oscillation probabilities. We 
discuss the resolution of the sign-degenerate solution by using analytical expressions for the 
probabilities and confirm the results by numerical calculations. In sec. [3] the method is illus- 
trated with simulations of realistic experimental configurations. As examples we combine the 
fi —>■ e super beam experiments SPL, T2HK, and NOz/A with a 7 = 100 beta beam (e —>■ fi), 
all operating in the regime of relatively small matter effect. A summary is presented in 
sec. H 
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2 Resolving the sgn(Am| 1 )-degeneracy with T-conjugated oscilla- 
tion channels 

In this section we discuss how the sign- degeneracy can be resolved with CP and T-conjugated 
channels by considering the location of the degenerate solution in the plane of 6*13 and Sep- 
Similar considerations can be found for example in refs. [14,15,19,20]. Let us depart from 
an expression for the appearance oscillation probability in the v e — z/ M sector, expanded to 
second order in the small parameters S13 = sin# 13 and a = Am^/IAm^l valid for constant 
matter density [21-23]: 



sin 2 A(l - ahA) 2 2 sin 2 A A 

'app « * «i3 « 23 Jl _ ahA y + <* Sln 2012 C 23 ^ 2 
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where L is the baseline, E is the neutrino energy, and V is the effective matter potential [1]. 
The signs a, t, h describe the effects of CP-conjugation, T-conjugation, and the neutrino mass 
hierarchy, respectively: 

The matter effect enters via the parameter A. It is clear from eq. (00) that in the case of large 
matter effect A > 1 the terms (1 — ahA) depend strongly on the type of the mass hierarchy, 
and for ah — 1 (neutrinos and NH, or anti- neutrinos and IH) a resonance is encountered for 
A=l[3). 

In the following we will focus on a different situation, namely the regime of small matter 
effect A <C 1. Numerically one finds for a matter density of 3 g/cm 3 

- U ' Uy VGeV J V2.5 x 10-3 eV 2 J " 1 J 

Furthermore, we concentrate on experiments operating at the first oscillation maximum, 
which is characterized by A ~ 7r/2, or 

E^0.2Gev(^—) ( |A ^il ) . (5 ) 

V 100 km; V2.5 x 10- 3 ev 2 y v ; 

Eqs. (Hj) and (jSJ) imply that for experiments at baselines of 130 km, 295 km, 730 km the 
matter effect is of order 2%, 5%, 13%, respectively0 Keeping this in mind it makes sense 



*A discussion of the issue in which regions of L and E the matter effect is important can be found for 
example in ref. [23] . 
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to expand the probability eq. (OQ) also in the small quantity A. To simplify the following 
equations we set 6*23 = vr/4 and use the abbreviation 

a = sin2 £ 12 _J^. (6) 

Then eq. (OQ) becomes to first order in A 

-P ap p ~ 2s^ 3 sin 2 A + - a 2 + 2h a si 3 sin A cos(/iA — at<5 CP ) 

+ 2a S13A (sin A — A cos A) [2h S13 sin A + d cos(/iA — ai^cp)] > (7) 

where the first line is just the vacuum probability and the second line corresponds to the 
leading order matter effect correction. 

The reason why it is difficult to determine the neutrino mass hierarchy is a parameter 
degeneracy with s 13 and 5qp [14,15], i.e., for given s 13 , <5 CP , and sgn(Am| 1 ) in many situations 
the same probability can be obtained for the opposite sign of Am^ and different values s' 13 
and S' CP : 

P app (a, t; h, S13, Sep) = P apP (a, t; -h, s' 13 , 5' CP ) . (8) 

Assuming a given oscillation channel t and neutrino and anti-neutrino data (a = ±1) this is 
a system of two equations for the variables s' 13 and 5' CP . If a solution to this system exists 
the mass hierarchy cannot be determined. For example, in the case of vacuum oscillations 
it follows immediately from eq. ([7D that the condition (jSJ) can be fulfilled for s' 13 = si 3 and 
$cp = n -$cp [14]: 

P™(a, t; h, sia, 6 CP ) = P™{a, t; -h, s 13 , vr - S CP ) , (9) 

where P™p i s given by the first line of eq. (J7|), and eq. © holds independently of a, t and 
the neutrino energy E. 

To include the leading order matter effect correction we introduce small deviations from 
this vacuum solution: 

s 'i3 = s i3(l + e«) , S' CP = n - 5 C p + e s , (10) 

with e s ,es 1. Using eqs. <^ and ( fTUj) in eq. (JHI) and expanding to first order in e s , es, and 
A yields the condition 

[e s — 2ahA(l — A cot A)] [2h Si 3 sin A + d cos(/iA — atScp)] = ate^a sin(/iA — at5cp) ■ (11) 

For a given "true" hierarchy h, a fixed oscillation channel t, and a = ±1 this is a linear system 
of two equations for e s and es, which in general has a unique solution. Hence, for neutrino plus 
anti-neutrino data in one oscillation channel the leading order matter effect cannot break the 
degeneracy. This is the reason why experiments at relatively small baselines (L < 700 km) 
have very poor sensitivity to the mass hierarchy^] In order to resolve the degeneracy one has 



2 Eq. (Ilip can be fulfilled exactly only for one energy. Hence, in principle spectral information can be 
used to resolve the degeneracy. Note, however, that this is an effect at third order in the small quantities 
e s , es, A, S13, a. Hence, it is difficult to obtain enough statistics to explore spectral information. 
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to enter the regime of large matter effect, where the non-linear character of eq. (CQ) becomes 
relevant and prevents a solution of the two conditions (jHJ) for fixed t and a ± 1. 

However, the immediate conclusion from eq. ( TTTT) is that if all CP and T-conjugate chan- 
nels are available one obtains four independent relations (corresponding to a = ±1 and 
t = ±1) for the two variables e s and eg. Obviously, in general such a system has no solution 
and hence the degeneracy is broken. To illustrate this explicitly let us consider for simplicity 
the case A = n/2, i.e., experiments exactly at the first oscillation maximum. Then eq. (TTTT) 
simplifies to 

(e s — 2ahA) (2si3 + ata sin Sep) = at es a cos Sep ■ (12) 
For a = ±1 and given t this system of equations has the solution 

, a . 
e s = ht A — sin Sqp , 

Sl3 



£5 = /it : Sin d C P 

cosdep \s 13 

The crucial observation from these expressions is that the signs of both, e s and e$, depend on 
the oscillation channel t. Hence, with increasing matter effect A the location of the solution 
with the wrong hierarchy in the #13 — Sep plane moves in opposite directions for \i — > e and 
e — > \i transitions. 

We have verified this behaviour by numerical calculations of the oscillation probability 
for constant matter (without any expansion in small quantities). Based on such calculations 
we graphically solve the system of equations (jSJ) in fig. [H The appearance probability is 
calculated for NH and the parameters sin 2 26* 13 = 0.02 and Sqp = 36° (marked as a star). 
Then the curves in fig. [1] correspond to the set of values sin 2 29[ 3 and S' CP leading to the 
same probability but for IH. In each panel there are four curves, corresponding to the four 
combinations of neutrinos/anti- neutrinos and (e — > /i)/(/i — > e). The dots show the location 
of the degenerate solutions, where the curves for neutrinos and anti-neutrinos for a given 
channel cross. The second crossing close to the original parameter values correspond to the 
so-called intrinsic degenerate solution [15,24] (with the wrong hierarchy). For experiments 
at the first oscillation maximum this degeneracy is resolved quite efficiently by spectral 
information, and hence, in many cases it does not appear as a viable solution (see, e.g., 
ref. [25] for an explicit discussion in the case of the T2HK experiment). Therefore we will 
neglect it in the present discussion and focus on the solutions marked with dots in fig. [TJ 

In the upper-left panel of fig. [1] we consider a (hypothetical) experiment at a very short 
baseline L = 10 km and an energy of 0.02 GeV, where the matter effect is negligible, see 
eq. 01]). One finds that all four curves meet in the point corresponding to eq. © with 
S' C p = 7r — Sqp, which makes the determination of the hierarchy impossible even if all CP 
and T-conjugated channels are available. By increasing the baseline (and simultaneously 
choosing the energy to stay in the first oscillation maximum) one observes from the plots 
that the degenerate solutions for the e — > \i and fi — > e channels separate and move in 
opposite directions, in agreement with eq. ffT3"]) . 
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Figure 1: Location of the sign-degenerate solutions for different baselines. The plots show a graphical 
representation of the equations Pn PP (a, t; NH, s%s, Sqp) — f a pp(a, t; IH, s' 13 , S' CP ) for a = ±1 (neutrinos/anti- 
neutrinos) and t = ±1 ((e — > — * e)). The star indicates the assumed values for sin 2 2#i3 and (5cp 
for the NH, whereas the axes correspond to the primed parameters for the wrong hierarchy. The dots show 
the location of the degenerate solutions, where the curves for neutrinos and anti-neutrinos for fixed t cross. 
The four panels correspond to different baselines, and the neutrino energy in each panel is determined by 
assuming the first oscillation maximum and jAmf^! = 2.5 x 10~ 3 eV 2 , see eq. ([5]). 

3 Numerical simulations 

In this section we show by numerical simulations of realistic experimental configurations 
how one can benefit from the synergy of T-conjugated oscillation channels. We consider the 
three super beam experiments SPL, T2HK, and NOz/A, providing the v ^ — ► v e information, 

as well as a 7 = 100 beta beam (BB) experiment operating in the v e — > v ^ channels. The 
simulation is performed with the GLoBES software [26], and all setups correspond to the 
pre-defined configurations provided by GLoBES 3.0. The most relevant parameters for each 
experiment are given in tab. [TJ BB and SPL are CERN based experiments, using the 500 kt 
water Cerenkov (WC) detector MEMPHYS at Frejus, at a distance of 130 km, details can be 
found in ref. [16]. With such a short baseline and low energies the matter effect is very small. 
T2HK is the second phase of the T2K experiment [27] in Japan, based on a 4 MW upgrade 
of the beam and the 500 kt HyperKamiokande detector. Further details of the simulation 
can be found also in ref. [25]. In this experiment the matter effect is somewhat larger than 
for the CERN-MEMPHYS configuration, but still too small to explore the neutrino mass 
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Table 1: Main parameters of the simulated setups [26]. For BB the "beam" column corresponds to the 
number of useful 18 Ne ( 6 He) decays per year, whereas for the super beams the beam power is given. The 
systematical error cr sys corresponds to the uncertainty on the signal and background rates, uncorrelated 
between signal, background, neutrinos, and anti-neutrinos. 



hierarchy. Finally, we consider the Fermilab based NOz/A experiment [28] with a 25 kt totally 
active scintillator detector (TASD) at a baseline of 812 km, where the matter effect starts 
to be important. Note however, that here we consider an initial stage setup for NOz/A, with 
significantly less statistics than the other configurations. 

To simulate the data we take the following values for the oscillation parameters: Am^ = 
7.9 x 10~ 5 eV 2 , sin 2 #12 = 0.3, Am^ = +2.4 x 10~ 3 eV 2 , sin 2 # 23 = 0.5, and we assume an 
external uncertainty of 4% on the solar parameters, and 5% on the matter density uncertainty 
along the baseline of each experiment. Note that each experiment includes neutrino and anti- 
neutrino data, and appearance and disappearance channels are used in the analysis. 

In fig. [2] we show an example of the sign-degenerate solutions in the sin 2 26^ — #cp plane 
for the four experiments of tab. [TJ This plot confirms the behaviour discussed in the previous 
section on probability level by performing an actual fit to simulated data: One can see that 
the best fit point with the wrong hierarchy moves in opposite directions for BB (e — > /i) and 
the super beam experiments (/i — > e), relative to the vacuum solution at 5qp = n — 5cp e , 
which is indicated by the dashed line. Furthermore, the dislocation of the degeneracy for the 
super beam experiments from the vacuum value of 5cp is proportional to the baseline, i.e., 
increasing in the order SPL, T2HK, NOzM.. This behaviour shows that the combination of 
fi — > e and e — > /i experiments offers a promising synergy in resolving the sign-degeneracy. 
Note that that in all cases only one solution appears with the wrong hierarchy, which justifies 
to neglect the intrinsic degeneracy in the discussion of the previous section. 

Let us stress that in each case a detailed quantitative study is necessary in order to fully 
assess the potential of this method. The final ability to disfavour the degenerate solution 
depends on many details which affect the size of the allowed regions for the individual 
experiments. For example, the relatively large allowed region for NOz/A seen in fig. [21 which 
is a consequence of the much smaller statistics compared to the other setups, will certainly 
limit the sensitivity, tough the best fit point of the degeneracy for NOz/A clearly follows the 
trend discussed above. For the SPL/BB combination the effect discussed here has apparently 
not be found in refs. [20,29] . It is clear from fig.[2]that for SPL+BB a slight enlargement of the 
allowed regions would be enough to corrupt the ability to resolve the sign- degeneracy. The 
size of the allowed regions depends rather sensitively on the details of the analysis (see e.g., 
fig. 6 of ref. [16] for the case of SPL). A possible reason for the differences with refs. [20,29] 
might be for example the assumptions on systematics (see tab.[T]), or the inclusion of spectral 
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Figure 2: Allowed regions at 90% CL in the sin 2 26 13 — Sep plane with the wrong hierarchy for BB, 
SPL, T2HK, and NO^A. Data is simulated for NH and sin 2 20f 3 ue = 0.03, 6^ c = 0.15tt, marked with a 
star. The dots show the best fit point with 1H hierarchy. The horizontal dashed line indicates the value 
Scp = T — #cp e , corresponding to the location of the degenerate solution in vacuum, see eq. (0. For the BB 
external accuracies on lAm^ and #23 of 5% and 10% have been assumed, respectively. 



information for appearance as well as disappearance channels, which has quite a significant 
impact on the size of the allowed regions [16]. 

In fig. [2] we have assumed a "true" NH. If the true hierarchy is inverted the degenerate 
solutions for the fi — > e and e — > \x channels move into the opposite directions as in the case 
of a true NH. This follows from eq. ( {TBI , where h corresponds to the true sign of Am^. Of 
course, the complementarity between the T-conjugated channels remains independently of 
the true hierarchy. 

Let us comment on the widely discussed strategy of using information from two // — > e 
experiments at similar L/E but different baselines, see, e.g., ref. [30]. Such a situation 
is realized by the combination of T2K and NOi^A [31] or by placing a second detector in 
the NOz/A beam-line at a suitable off-axis angle, as proposed in ref. [32]. It is clear from 
the preceding discussion and from fig. [2] that the synergy from such a combination is less 
effective than using e — > fi information. The reason is that the degeneracies for two // — > e 
experiments move in the same direction in the sin 2 2# 13 — <5cp space, only the size of the 
dislocation is different due to the different baselinesjf] 

3 Note that the present discussion does not necessarily apply to super beam experiments with neutrino 
energies E ~ GeV and a detector at large baselines L > 1000 km, far beyond the first oscillation maxi- 
mum [11, 12]. Apart from the fact that in this case the matter effect cannot be considered to be small, also 
the rich information from the energy spectrum at higher oscillation maxima can lead to a rather different 
behaviour of the sign-degenerate solution than discussed here. 
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Figure 3: Minimal value of sin 2 2$i 3 for which the IH can be excluded at 2a (A% 2 = 4) if the true 
hierarchy is normal, as a function of Sep (left) and the fraction of all possible values of <5cp (right). The 
dashed lines correspond to the super beam experiments SPL, T2HK, NO;/A, whereas for the solid lines each 
of these super beams is combined with the beta beam. For comparison also the combination NOz/A+T2HK 
is shown. 

In fig. [3] we show how the combination of the super beam experiments SPL, T2HK, 
and NOz/A with the 7 = 100 beta beam significantly enhances the sensitivity to the mass 
hierarchy due to the (// — > e)/(e — > fi) synergy. The dashed curves show the sensitivities of 
the super beams alone. One observes that these experiments can assess the mass hierarchy 
only in a certain range of 5qp values, and there is no sensitivity even for sin 2 2#i3 = 0.1 for 
75% (50%) of all values of 5 CP for SPL (T2HK, NOz/A). Note that NOz/A and T2HK have a 
rather similar sensitivity, where NOz/A has the advantage of the longer baseline of 812 km, 
whereas the short baseline of T2HK of 295 km is compensated by the large statistics implied 
by the 4 MW beam and 500 kt detector. The main limitation of SPL is of course the short 
baseline of 130 km. 

When these experiments are combined with the BB the sensitivity is significantly im- 
proved, see solid lines in fig. [SJ The dependence on 5np practically disappears and a sta- 
ble sensitivity is obtained for sin 2 2#i3 > 0.02 — 0.03jj The effect is most remarkable for 
SPL+BB [16], since none of these experiments on its own has any notable sensitivity. Indeed, 
in the parameter range shown in fig. [3] there is no sensitivity for BB alone. The sensitiv- 
ity of the SPL+BB combination fully emerges from the complementarity of T-conjugated 
channels in the small matter effect regime. Also for T2HK and NOz/A the effect is clearly 
visible, and again the larger baseline of NOz/A is compensated by statistics in T2HK. For 
comparison we show also the combination NOz/A+T2HK (without BB) in fig. [3] (see ref. [31] 
for detailed discussions of the case NOz/A + T2K phase I). Also in this case the sensitivity 

In fig. [3] we have assumed that the true hierarchy is normal, but we have checked that for a true IH the 
results are very similar. 
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improves, however the complementarity is much less pronounced, a dependence on Sqp re- 
mains, and the effect is more similar to the addition of statistics than a true synergy (see 
also the corresponding discussion related to fig. [2]). 

As a side remark let us mention the possibility pointed out in the second paper of ref. [7], 
that very accurate measurements of the neutrino mass-squared difference in v e and disap- 
pearance experiments allow in principle to distinguish between NH and IH (even in vacuum). 
In the beta beam/super beam combination both disappearance probabilities P ee and P w are 
observed. However, we have checked that for the experiments considered here numerically 
this effect is completely negligible, and practically the full sensitivity shown in fig. [3] emerges 
from the matter effect in the appearance channels. 

4 Summary and discussion 

In this letter we have considered neutrino oscillation experiments operating at the first 
oscillation maximum in the regime of small matter effect, i.e., at relatively short baselines 
of several hundred km. In such a case there is very poor sensitivity to the neutrino mass 
hierarchy, if only data from neutrinos and anti-neutrinos are available in a fixed oscillation 
channel. The reason is that the leading order correction in the small matter effect parameter 
A cannot break the sgn(Am| 1 )-degeneracy. The usual strategy to resolve the degeneracy is 
to enter the regime of strong matter effect (by going to longer baselines and higher neutrino 
energies), where the non-linear dependence of A becomes important. Here we have proposed 
an alternative method, namely the combination of all four CP and T-conjugated oscillation 
channels. We have shown that the location of the sign-degenerate solution in the plane 
of #i3 and <5cp moves in opposite directions for the v e and V e — »V M channels when 
the matter effect increases. This synergy allows to resolve the degeneracy even for matter 
effects as small as a few percent. We have discussed the method at the level of oscillation 
probabilities, and illustrated the effect also by simulations of representative experimental 
configurations. 

A typical situation where our method applies is the combination of beta beam (e — > /i) 
and super beam (// — > e) experiments. We have demonstrated that a significant synergy 
exists for the determination of the mass hierarchy by simulations of the SPL, T2HK, and 
NOi^A super beams combined with a 7 = 100 beta beam. Most remarkable, even for the 
SPL /beta beam combination, where both experiments have a baseline of only 130 km, there 
is sensitivity to the mass hierarchy for sin 2 2#i3 > 0.03 due to the synergy of the T-conjugated 
channels. Another possibility to take advantage of this effect could be a low energy neutrino 
factory [17] operating in an L/E regime where the matter effect is not yet fully developed. 
The liquid argon technology, which has been considered for the detector in that reference, 
has excellent sensitivity for muon as well as electron detection. If in both cases the charge 
can be identified all four oscillation channels and V e — >V M were available at the 

neutrino factory. If charge identification is not possible for electrons one could combine the 
neutrino factory with a super beam experiment providing the fi — > e information. 
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The purpose of this note is to point out the existence of a synergy of T-conjugated chan- 
nels, which can significantly increase the sensitivity to the neutrino mass hierarchy. Whether 
this method for the mass hierarchy determination is indeed useful for a given experimen- 
tal configuration, or is competitive with alternative approaches needs to be confirmed by 
detailed simulations and comparison studies, which is beyond the scope of this work. 
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